Dynamical properties of the one-dimensional spin- 1/2 Bose-Hubbard model near 
Mott-insulator to ferromagnetic liquid transition 
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We investigate the dynamics of the one-dimensional strongly repulsive spin-1/2 Bose-Hubbard 
model for filling v < 1. While at v — 1 the system is a Hubbard-Mott insulator exhibiting dynamical 
properties of the Heisenberg ferromagnet, at u < 1 it is a ferromagnetic liquid with complex spin 
dynamics. We find that close to the insulator-liquid transition the system admits for a complete 
separation of spin and density degrees of freedom valid at all energy and momentum scales within 
the t — J approximation. This allows us to derive the propagator of transverse spin waves and the 
shape of the magnon peak in the dynamic spin structure factor. 



Collective dynamics of interacting quantum particles 
in a periodic potential is a paradigmatic problem of 
many-body physics. The simplest model incorporating 
the competition between the interaction and kinetic en- 
ergy on a lattice is the Hubbard model [H . For fermions, 
this model has been the cornerstone of the theory of 
strongly correlated systems. In particular, it describes 
the transition from a Mott insulator at a commensurate 
filling to a metallic state when the system is doped. The 
nature and properties of the weakly doped state remain 
the subject of fierce debate, in which an interplay be- 
tween the spin and charge dynamics occupies a central 
place. It is this interplay that is believed to lead to a rich 
phase diagram with putative phases ranging from a "nor- 
mal" Fermi liquid state with heavy holes propagating in 
an antiferromagnetic background to a d-wave supercon- 
ducting state arising from a spin liquid. The bosonic 
version of the Hubbard model, the Bose-Hubbard (BH) 
model 0, also has a wide range of applications, from 
granular superconductors to Josephson junction arrays. 
Recent experimental progress in ultracold atomic gases 
further increased the interest to this model 0]. 

Like its fcrmionic sibling, the BH model exhibits a 
Mott-insulator to liquid transition. This transition has 
already been observed experimentally in the ultracold 
one-component ("spinless") boson gases 0, 0, Q and 
is likely to be seen in multicomponent ("spinful") gases 
soon. A principal difference between the fermionic and 
bosonic case comes from their different magnetic orders: 
while the fermionic ground state is antiferromagnetic, 
the bosons have a tendency to ferromagnetism. This 
difference should have important consequences for the 
Mott-insulator to liquid transition. The essentially non- 
perturbative nature of this transition makes its theo- 
retical analysis a non-trivial task. Therefore, the par- 
ticular case of one special dimension (ID), where non- 
perturbative methods are well developed is of special im- 
portance. For both fermionic Hubbard and spinless BH 
models in ID the phase obtained by doping the Mott 
insulator is the Luttinger Luquid (LL), whose properties 



are well established 0. However, due to ferromagnetism, 
the doped state of the BH model with spin cannot be a 
LL. For a continuous system it was demonstrated very 
recently Q that ferromagnetism leads to a new phase, 
named ferromagnetic liquid, with very unusual proper- 
ties 0, 0, [HI [HI ■ What is the analogue of that phase 
for a doped BH insulator, particularly in the vicinity of 
the insulator to liquid transition, is an open question. 
The analysis of this problem is especially deserving since 
spin-resolved measurements in a ID boson quantum gas 
were already reported [ljj . 

In this Letter we consider the two-component (spin- 
1/2) ID Bose-Hubbard model with strong on-site re- 
pulsion. We find that in the vicinity of the liquid to 
Mott-insulator transition the system admits for a com- 
plete spin-charge separation at all energy and momen- 
tum scales within the t — J approximation. This result 
makes possible explicit calculation of dynamic correlation 
functions of the system. In particular, we calculate the 
Green's function of transverse spin excitation (magnon) 
and derive the shape of the magnon peak in the dynamic 
spin structure factor. 

The Hamiltonian of the repulsive (U > 0) spin 1/2 BH 
model on a ID lattice with M sites is 
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, b a j and Qj 



bijbij+b^jb^j are the 



In this expressions b a ^ 

creation, annihilation and the particle number operators 
on the j-th lattice site, respectively. For a given particle 
number N one defines the filling factor v = N/M. 

We derive an effective Hamiltonian valid for large on- 
site repulsion, U ^> th, from Eq. {T]) in analogy with the 
t-J approximation for the Hubbard model [1 31 ] - This is 
done by treating the potential term in Eq. (JT|) as the 
leading order Hamiltonian, the kinetic term as a per- 
turbation, and expanding over this perturbation to the 
second order. It is convenient to write the resulting t-J 
Hamiltonian using the operators c\ c and £, which have 
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been previously used for the analysis of the collective dy- 
namics of spin and charge excitations in fermion [l4| and 
boson [1, [13] systems. We have H e g = H c + H s , where 



H c = -th X]( c i c J+ 1 + h - c -)' 
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and 
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t +lCj t Cj+1 c 3 -) x [^-) + + l)] 2 . (3) 



In these expressions (cj) are spinless fermion creation 

(annihilation) operators, = ejej, and A/} = X)n=i & n 
counts the number of particles between the first and j-th 
lattice site. The operator £ s (Afj), s = x, y, z represents 
the spin state of the nearest particle to the left of the site 
j and is defined by 
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where £ s (m) are local spin-1/2 operators on a lattice, 
[£ x (m),£ y (n)] — i5 mn £ z (n). Note that in the sector with 
maximal value of the total spin Eq. ([3]) gives the t — J 
approximation of the spinless BH model, cf. Eq. (4) of 
Ref. [H. 

While H c generates the time evolution of free spinless 
fermions, spin and charge (density) degrees of freedom 
are coupled in H s non-trivially. An advantage of the rep- 
resentation (J3j) is that in the vicinity of the Mott phase, 
1 — v <C 1, which we focus on in our Letter, one can ne- 
glect quantum fluctuations of c-dependent operators in 
H s . As a result one gets the Hamiltonian of the Hciscn- 
bcrg fcrromagnct 
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H s = -2JJ2 



WU + 1) 



(5) 



with the exchange constant J = t"i/U, We shall see that 
already within the approximation (J5]) the dynamics of 
the system is non-trivial. The analysis of perturbative 
corrections in 1 — v to Eq. ([5]) remains an interesting 
open question. 

We now proceed to the investigation of spin dynamics 
above the fully polarized ground state | ft) of the Hamil- 
tonian (fTJ. Like the ferromagnetic Bose gas H, a non- 
trivial part of this dynamics is encoded in the correlation 
function of transverse spin excitations 



Gi(j,«) = a|s + (j,t)s-(0,0)|t), 



(0) 



where s + (j) = £>L-bjj and s_(i) = [s+(j)]L Since 
G ± (j,t) = G±(-j,t) and G* ± {j,t) = G±(-j, -t) we con- 
sider j, t > 0. We focus on the regime of weak hole doping, 



/Lt, and strong on-site repulsion: 

fj, = 1 - v < 1 and Ji,; 1 = t h U^ x < 1. (7) 

Length and time scales in our problem in the regime ([7|) 
are set up by the parameters qf and tp, respectively: 
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Ef — <lpth- 



(8) 



Combining Eqs. ([7]) and ||SJ) we get three dimensionless 
small parameters: qF, J/th, and Ep/tt- 

To calculate G± using Eqs. and ((5]), we express the 
operators s± through c\ c, and t± = £ x ± i£ y : 



(9) 



and factorize the ground state | ft) = | •ff)j? ® |FS), where 
ft) a is the fully polarized ground state of H s , Eq. (O, 
and |FS) is the ground state of H c , Eq. (|2|), at a given v. 
For Eq. ([6|) this gives 



G±(j,t) = / d\G H {\t)D v {\- 3 ,t). 



(10) 



Here Gn(^,t) = exp{— iE s {\)} /2-k is the space Fourier 
transform of the transverse spin-spin Green's function 
H (ft \£+(n, t)£-(0, 0)| ft) of the Heisenberg ferromag- 
nct ([5|), and E S (X) = 2Jt(l — cos A) is the dispersion of a 
single magnon in that model. The fluctuations of density 
enter Eq. (p~0|) through the Green's function [l6| 



D v (\;j, t) = (FSIe-^We^^lFS). (11) 

We stress that though Gh {D v ) contains collective spin 
(density) degrees of freedom only, spin and density arc 
coupled to each other in G± through the integration over 
A in Eq. ([TO]). 

There are two ways to analyze Eq. (fTTj) . Like other free 
fermion Green's function s 11711 . it can be represented as 
a Fredholm determinant [18| |. The large space, qpj > 1, 
or time, t/tp > 1, asymptotics of this determinant can 
be found using the results of Ref. [l9| , though the calcu- 
lations are rather involved. Alternatively, large space- 
time asymptotics of Eq. (fTTj) can be calculated using 
the LL formalism. There the operator Afj(t) is given 
by Nj(t) = vj — n~ 1 d x (j)(x,t), where is a free bosonic 
field and x = aj, a being the lattice spacing (2fJ]. Thus 



D„(\;j,t) ~ exp 



J F"c 



Here vf = 2thqF is the Fermi velocity of the holes, which 
coincides with the sound velocity of the system. The cut- 
off parameter t c ~ <f ■ Its explicit value cannot be deter- 
mined within the LL formalism, while the analysis of the 
Fredholm determinant gives t c = e _1_7_ln4 tF ~ 5.16 x 
10~ 2 £f, where 7 = 0.577 ... is the Euler-Mascheroni con- 
stant 2l|. Having the result (Tj"2"|) . we can investigate 



different parametric regimes shown in Fig. [TJ 
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FIG. 1: Shown are the parametric regimes (A), (B), and (C) 
corresponding to different relative values of the parameters 
qf — Tfi, J/th, and Ep/th- 



(A): J/th > qp. The doping is so small that during 
a large time t* the spin excitation evolves like in the 
undopcd system: 



G ± (j,t)=e- 



t < t* 



(13) 



where Jj is the Bessel function of the first kind. To 
estimate t* note that Eq. (fT3|) can be obtained by letting 
v = 1 in Eq. dTTJ) which implies D v= i (A; j, t) = e iX] . Such 
an approximation is valid as long as qpj < 1 and t < tp. 



As j —> oo one has Jj ~ e 



-3 1" 3 



starting from j ~ 2 Jt and 



due to the exponential decay possible deviations of G± 
from what is predicted by Eq. (| 13[) become immaterial 
for larger j. One thus arrives at the condition 2Jqpt < 1, 
which implies i» = 1/2 Jqp < tp. 

For t > tp one can use Eq. (fT2|) and the integral 
in Eq. (fT0|) is dominated by the saddle point A s = 
arcsin (— j/2 Jt — iO). Therefore 



i Jt cos A, 



■G H (X s ,t)D v (X s ;j,t) (14) 



with D v given by Eq. (fT2"|) . One can see that G± oscillates 
as a function of j for j < 2Jt and decays as e - -' 1113 for 
j > 2 Jt. The divergence of Eq. (|14[) as j — » upi, custom- 
arily called light-cone singularity, indicates what fraction 
of a collective excitation propagates with the velocity of 
sound. 

(B) : Ep/th < J/th < qp- Here the result flTSJl is valid 
for t < tp and it crosses over to Eq. p4[) for t > tp. 
The main difference from the case (A) is the absence of 
the light-cone singularity. Indeed, the exponential de- 
cay of Eq. (fl~4"|) due to the imaginary part of A s starts 
before j reaches the light cone, Vpt. Note that one can 
use D v (X s ;j,t) ~ e~ iX ° j D V {X S ] 0, t) in Eq. dHJ) in this 
regime. 

(C) : J/th < Ep/th- Here for t < tp the spin excita- 
tion does not propagate at all, G±(j, t) = 8oj. For t > tp 
and 27T 2 Jt < \n(t/t c ) (the latter inequality could be ap- 
proximately written as t < J^ 1 ) the integral in Eq. (jTUJ) 
is dominated by the region |A| < 2tt/ y/ln(t/t c ) and the 
result of the integration is [22j 



Gx(j,t) 



21n(t/t c 



■ cxp 



2\n(t/t c ) 



(15) 



For this calculation we replaced \j — vpt\ by vpt in 
Eq. (fl"2"|) : keeping the light cone singularity in Eq. (fTS)) 
would be unphysical since Gx decays as a gaussian at 



distances j ~ y/]nt/t c <C vf^- The physical interpreta- 
tion of Eq. (fT5|) is that the propagation of spin excitation 
is caused by the fluctuations of the position of the parti- 
cle carrying this spin rather than by exchange processes. 
The same propagation mechanism takes place in the con- 
tinuum limit in the trapped regime, cf. Eq. (6) of Ref. 
Finally, for t > tp and 2n 2 Jt > ln(t/t c ) the asymptotics 
of G± is given by Eq. (|14|) and the analysis carried out 

as 



for the case (B) applies. In particular, G± ~ e - ^ 3 
j — > oo, which is pronouncedly different from the gaus- 
sian decay of Gj_ in the continuum limit, cf. Eq. (16) 
of Ref. Q . This difference stems from different propaga- 
tion mechanisms for the spin excitations in the v< 1 and 
1 — v <C 1 regimes: while in the former both longitudinal 
and transverse spin excitations are carried by particles, 
in the latter longitudinal spin (that is, density) excita- 
tions are carried by holes and transverse spin excitations 
by particles. 

The rest of this Letter is devoted to the anal- 
ysis of the spectral function A{k,w) = ^2^8{lo — 
E f (k))\(f,k\s.(k)\ ft)| 2 , where H\f,k) = E f (k)\f,k) 
and / enumerates all the states with a given momen- 
tum k. Using Eq. pTT)) and the identity D v (X;j,t) = 
e _lAj Z) ft (— A; j, t) following from particle-hole symmetry 
we write 

A{k,uj) = J ^D„(X;k-X,u-E s (X)), (16) 

where is the Fourier transform of D^(X]j, t) with re- 
spect to j and t and E S (X) is defined below Eq. (fT0|) . For 
k and w — * Eq. (fTS|) implies 

Dfj,(X; k,u) ~ e-^9(Lo - 2tt\ P \)[u 2 - (2np) 2 ] a ^ , (17) 

wherep = k + A/i and a(X) = — l + (A/27r) 2 . Substituting 
Eq. (fl~7l) into (p~6|) we arrive at 



A(k,u) ~ 0(u - E s {k/v))[uj - E s {k/v)]^ k) (18) 

valid for ui close to E s (k/v) and — it < \k\ < it. The 
threshold exponent A(fc) in Eq. ([18]) is 



(19) 



Equations (|18[) and (flU)) constitute one of the main re- 
sults of this Letter, so let us discuss them in more detail. 
The threshold exponents of the spectral functions were 
calculated for several classes of ID interacting models, 
mostly using quite involved techniques (23l |. Probably, 
the closest to our work are the studies of A(k,u) in ID 
ferromagnetic Bose gas, the model defined by taking the 
continuum limit of Eq. ([1]). In this model A(fc), found in 



Refs. [alia L!J|, is also given by Eq. ([19]), but with differ- 



ent coefficient in front of A; 2 term. Moreover, according to 
Ref. [ll[ , any Galilean-invariant system of Bose particles 
with ferromagnetic ground state should have A(fc) given 
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FIG. 2: Shown is the spectral function A(k,u>) as a func- 
tion of ui for several values of k: 3tt/8 (red), 5n/8 (green), 
37r/4 (blue), and ir (pink). The curves are obtained from the 
numerical analysis of the Fredholm determinant in Eq. (|16|l . 
The positions of the singularities following from the asymp- 
totic expression (|18[) are indicated with vertical dotted lines. 
The parameters of the model used for getting the plot are 
Jt,; 1 = 0.04 and q F = 0.1, which implies E F /2J = 0.125. 



by Eq. (jT9|) . We see that the BH model, which is not 
Galilean- invariant, still has the threshold exponent (fT9]l 
at 1 — v <C 1. Whether this is an accident, or a manifes- 
tation of some fundamental principle other than Galilean 
invariance, is an open question. A deeper understanding 
of this issue could be obtained by calculating the thresh- 
old exponent for the BH model at arbitrary filling, which 
is not yet done. 

To go beyond the asymptotic result (fTS)) we per- 
form [211 ] the numerical analysis of the Fredholm determi- 
nant representation of Eq. (fTTj) , given in Ref. [HI . The re- 
sulting plots of A(k,ui) for Jt^ 1 = 0.04 and qF = 0.1 arc 
given in Fig. [2l One can clearly see in this figure that (i) 
The position of the threshold singularity is described by 
Eq. (fT8|) . (ii) A(k,u) extends below the threshold, which 
is particularly well seen for k = 3ir/4. This effect is due to 
the emission by a spin excitation of multiple particle-hole 
pairs with the dispersion lying below E s (k/v). Within the 
numerical precision we observe an abrupt cutoff of the 
spectral weight at some finite frequency below E s {k/v) 
(for example, at oj/2J = 1.6 for k = 37r/4). The rea- 
son for such a behavior is unclear, (iii) To the right of 
the threshold singularity a second peak develops with 
increasing fc, becoming very pronounced for k — it. The 
theoretical explanation of its appearance is an open prob- 
lem. 

In conclusion, we have investigated the dynamics of 
the spin-1/2 Bose- Hubbard model with strong on-site re- 
pulsion and small hole doping. We showed that this sys- 
tem exhibits a transition from a Mott-insulator phase to 
a ferromagnetic liquid. We found collective variables in 
which spin and charge dynamics separate at all energy 
and momentum scales within the t — J approximation. 
Using these variables we obtain the transverse spin-spin 



correlation function G±, Eq. ([6]), in the form (fT0|) and 
(jlip . convenient for both analytic and numerical investi- 
gation. For the spectral function A(k, oj) associated with 
G±, we have found the analytic expression for the shape 
of the magnon peak, Eqs. (fig)) and (fTi?)) . Numerical anal- 
ysis of Eq. (fTfj)) revealed a non-trivial fine structure of 
the magnon peak, as shown in Fig. [2] 
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